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Abstract 

We study the two-dimensional Wess-Zumino model with extended 
N = 2 supersymmetry on the lattice. The lattice prescription we 
choose has the merit of preserving exactly a single supersymmetric in- 
variance at finite lattice spacing a. Furthermore, we construct three 
other transformations of the lattice fields under which the variation of 
the lattice action vanishes to 0{ga?) where g is a typical interaction 
coupling. These four transformations correspond to the two Majorana 
supercharges of the continuum theory. We also derive lattice Ward 
identities corresponding to these exact and approximate symmetries. 
We use dynamical fermion simulations to check the equality of the 
massgaps in the boson and fermion sectors and to check the lattice 
Ward identities. At least for weak coupling we see no problems asso- 
ciated with a lack of reflection positivity in the lattice action and find 
good agreement with theory. At strong coupling we provide evidence 
that problems associated with a lack of reflection positivity are evaded 
for small enough lattice spacing. 



1 Introduction 

Supersymmetry is thought to be an important ingredient of many theories 
which attempt to unify the separate interactions contained in the standard 
model of particle physics. Since low energy physics is manifestly not su- 
persymmetric it is necessary that this symmetry be broken at some energy 
scale. A set of non-renormalization theorems ensure that if SUSY is not 
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broken at tree level then it cannot be broken in any finite order of pertur- 
bation theory see eg. Q. Thus we are led to investigate non-perturbative 
mechanisms for SUSY breaking. The lattice furnishes the only tool for a 
systematic investigation of non-perturbative effects in field theories and so 
significant effort has gone into formulating SUSY theories on the lattice . 

Unfortunately, there are several barriers to such lattice formulations. 
Firstly, supersymmetry is a spacetime symmetry which is generically bro- 
ken by the discretization procedure. In this it resembles Poincare invariance 
which is also not preserved in a lattice theory. However, unlike Poincare 
invariance there is usually no SUSY analog of the discrete translation and 
cubic rotation groups which are left unbroken on the lattice. In the lat- 
ter case the existence of these remaining discrete symmetries is sufficient to 
prohibit the appearance of relevant operators in the long wavelength lattice 
effective action which violate the full symmetry group. This ensures that 
Poincare invariance is achieved automatically without fine tuning in the con- 
tinuum limit. Since generic latticizations of supersymmetric theories do not 
have this property their effective actions typically contain relevant super- 
symmetry breaking interactions. To achieve a supersymmetric continuum 
limit then requires fine tuning the bare lattice couplings of all these SUSY 
violating terms - typically a very difficult proposition. 

Secondly, supersymmetric theories necessarily involve fermionic fields 
which suffer from so-called doubling problems when we attempt to define 
them on the lattice. The presence of extra fermionic modes furnishes yet an- 
other source of supersymmetry breaking since typically they are not paired 
with corresponding bosonic states. Furthermore, most methods of eliminat- 
ing the extra fermionic modes serve to break supersymmetry also. 

In this paper we employ a lattice formulation of the two-dimensional 
Wess-Zumino model which was first written down in [^, In these earlier 
works the lattice formulation is found by discretizing the Nicolai map for the 
model . In our case we rederive the formulation in a slightly different way 
- we start from a simple discrete model which exhibits a one parameter local, 
supersymmetric invariance and show how this model may be generalized to a 
two dimensional Euclidean lattice field theory provided certain integrability 
conditions are satisfied. The N = 2 Wess-Zumino model is then found as the 
essentially unique solution to these conditions. Since the continuum model 
contains two Majorana supercharges we would expect the lattice model to 
possess three further transformations which are invariances of the action in 
the naive continuum limit. We construct these transformations explicitly 
and from them derive a set of exact and broken lattice Ward identities. 

To check these ideas explicitly we have simulated the simplest realization 
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of the model for a range of masses and couplings, computing both boson 
and fermion massgaps and the correlation functions needed for checking 
the supersymmetric Ward identities. To perform the simulations we have 
replaced the fermionic fields by commuting pseudofermionic fields in the 
usual manner. 

The outline of the paper is as follows; first we introduce a simple discrete 
model with an exact SUSY-like symmetry, showing how it can be used to 
describe a lattice version of supersymmetric quantum mechanics and then 
discussing its extension to two-dimensional field theory. The Ward identities 
are then introduced and we show how the expectation value of the total 
action (including the contribution of pseudofermion fields) can be used as 
an order parameter for SUSY breaking. Following from this theoretical 
introduction we present our numerical results both for weak and strong 
coupling. The final section contains our conclusions. 

2 Simple SUSY Lattice Model 

Consider a set of P real commuting variables Xi and two sets of P real grass- 
mann variables '0j and V'j with i = 1 . . .P governed by an action S{x,ilj,'ip) 
of the form 

S=^Ni{x)Ni{x)+lp,^^Pj (1) 

with the field Ni{x) an arbitrary function of Xj. It is easy to see that this 
action is invariant under the following SUSY transformation. 

SiXi = ViC 

SiA = 



5S = N.Q^^x,+5^,—^, 

which vanishes on account of the grassmann nature of the infinitesimal pa- 
rameter ^. Notice that the variation of the matrix 
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also vanishes for similar reasons. 

Let us now choose the fields x^ijj^ijj to lie on a spatial lattice equipped 
with periodic boundary conditions and take the fermion matrix Mij = 
to be of the form 

M,,=Df^+Pl'^{x) 

The symmetric difference operator Df^ replaces the continuum derivative 
and can be written in terms of the usual forward and backward difference 
operators. 

and Pij{x) is some (local) interaction matrix polynomial in the scalar fields 
X. The resulting model is easily recognized as supersymmetric quantum 
mechanics regularized as a + 1 dimensional Euclidean lattice theory [^. 
Furthermore it is trivial to find a field Ni [x) which yields this fermion matrix 
under differentiation 

iV, = Df^x, + Pl{x) 

Notice, however, that the resulting bosonic action \Nf is not a simple dis- 
cretization of its continuum counterpart 

^font = / dT\ [{drxf + {P' {x)f) 

as it contains a new cross term C = P'i{x)DfjXj which would be a total 
derivative (and hence zero) in the continuum but is non-vanishing on the 
lattice and required to ensure the transformation eqn. ^ is an exact symmetry 
of the theory. Notice that this extra term also vanishes on the lattice for a 
free theory where P'{x) = mx because of the antisymmetry of the matrix 

Notice that if I imagine changing variables in the partition function Z = 
J Dxe~^^^''^'^') from x to the Jacobian resulting from this transformation 
cancels the fermion determinant yielding a trivial gaussian theory in the 
field A^. This is an example of a Nicolai map and the existence of such a 
transformation of the bosonic degrees of freedom can be shown to imply an 
exact supersymmetry |^]. While most supersymmetric theories admit such 
a map, in the generic case it is non-local - that is the mapped Nicolai field 
N will be a function of arbitrarily high derivatives of the original boson field 
X. In the case of SUSY quantum mechanics (and as we will see later the 
N = 2 Wess-Zumino model) the expression is local. It can then serve as a 
basis for constructing a lattice theory with an exact supersymmetry as was 
pointed out in g, [| and |]. 
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So far we have neglected the fact that the form of the fermion action 
appears to admit doubles - the symmetric difference operator Dfj behaves 
like sin/ca in lattice momentum space yielding zeroes at both ka = and 
the Brillouin zone boundary ka = vr. Indeed, both the fermionic and bosonic 
actions now contain spurious modes which are not part of the continuum 
theory. The extra bosonic modes arise from using D^D^ as the kinetic op- 
erator rather than the usual scalar lattice Laplacian □ = Dj^D_. However, 
we can use our freedom in choosing the interaction matrix Plj{x) to add a 
Wilson term to the fermion action 

Plj{x) = —D^ + local interaction terms 

where the matrix = h(D'^- — Dj'-) = Dj,-. By construction this eliminates 
the doubles from the free fermion action completely; what is, perhaps, more 
surprising is that it also renders the boson spectrum double free too. This 
can be seen to be a consequence of the lattice supersymmetry. 

One further observation is in order. Consider a second supersymmetry 
transformation 

S2Xi = ViC 
62i^i = 



where 



dxi 



= -M^ = - Pl!^{x) 



The action in eqn. || is no longer invariant under this transformation 

262C (2) 



628 = -52{N?-N^ 



but transforms into the supersymmetry variation of (twice) the cross term 
C. As we have argued, for a free lattice theory or in the naive continuum 
limit this term will vanish and the model will be invariant under this sec- 
ond supersymmetry. For the lattice theory in the presence of interactions 
{P'{x) ~ gx'^,n > 1), this second symmetry will be broken by terms 0{ga'^) 
where the suppression by two powers of the lattice spacing reflects the fact 
that = dcont + 0(a^). Thus the second supersymmetry is broken only by 
irrelevant operators. Since quantum mechanics is a finite theory we then ex- 
pect that the continuum theory will have the two invariances that we expect 
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of supersymmetric quantum mechanics Q. We have verified this expUcitly 
in Q in which a computation of both the mass spectrum and the supersym- 
metric Ward identities revealed the existence of = 2 supersymmetry in 
the continuum hmit. 



3 The Lattice Wess-Zumino Model 

The action eqn. |l| and supersymmetry transformations eqn. |2| do not depend 
strongly on the existence of a background lattice of given dimensionality - 
indeed this physical interpretation only arises when we choose the form of 
the fermion operator. This allows us to use it as a basis for constructing 
candidate lattice field theories in higher dimensions which admit supersym- 
metry. 

In two dimensions the fermions will be represented by two independent 
two-component spinors whose components we will assume to be real (this 
restriction will turn out to be valid for N = 2 theories in Euclidean space). 
Thus we will imagine that the indices i,j can be promoted to compound 
indices i ^ i, a, j ^ j, j3 labeling spacetime and spinor components respec- 
tively. We immediately realize that there will be two scalar fields now in 
the theory Xi xf and the fermion matrix will take the form Dij — > D'^^ 
(from now on we will use D in place of D^). To maintain contact with the 
simple, discrete model we will require a Euclidean fermion operator which 
is also entirely real. Then the most general fermion matrix respecting this 
condition takes the form 

where A{x) and B{x) are real matrix fields and we have chosen a Majorana 
basis for the Dirac matrices so that 7^,7^ and ^7^ are also real. 
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72 = I 











71 = 

To remove the doubles we again add a Wilson term to the interaction matrix 
Aij = —Df- + interactions where 

The resultant fermion matrix is easily recognized as a discrete version 
of the continuum Wess-Zumino model and is the same fermion operator 
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appearing in |^],@|. Having chosen this fermion matrix we can attempt to 
find a vector A^" whose derivative yields M^j^ . Clearly must have the 
form 

where /° which represent mass and interaction terms must still be deter- 
mined. Ignoring for a moment the spacetime indices it is clear that strong 
restrictions are placed on the vector f". We must have 

^ ^ 9/^ ^ df ^ ^ 9/2 ^ _dfl 
dx^ dx'^ dx^ dx2 

Of course these are just Cauchy-Riemann conditions. In other words the 
integrability condition that M be a derivative of some vector imposes a 
complex structure on the scalar fields in the theory. Indeed, the bosonic 
part of the action can now be rewritten in terms of a complex vector t]^^^ {(p) 
whose real and imaginary parts are just the two components and A^^ 
respectively (we have again suppressed spacetime indices for clarity) Sb = 
ifjWfjW where Rec/) = and Im0 = x'^ and 

r^W = D,'^ + W'{cf)) 

where we have introduced complex coordinates z = {x+iy)/2, z = {x — iy)/2 
so that 

D, = Di- iD2 

with Di, Z?2 derivative operators in the two dimensional lattice. The sig- 
nificance of the superscript on rj^^^ will be become apparent later. W'{(p) 
is an arbitrary analytic function of the complex field with its complex 
conjugate. Furthermore, in this language the fields A and B are nothing 
but the real and imaginary parts of VF" ((/)). Expanding the bosonic action 
yields 

2,2 

The first two terms go over as a ^ to the bosonic part of the continuum 
action for the N = 2 Wess-Zumino model while the last two terms are 
clearly total derivatives which will vanish both in the continuum and for a 
free lattice theory. For an interacting theory they are necessary to preserve 
the lattice supersymmetry transformation. However they spoil the reflection 
positivity of the lattice action, a point we shall return to when we present 
our numerical results. 
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So far we have shown that the lattice action 

S=\NtNt+TiM^f^l 



admits the following invariance 



<5i< = i>ti 

Sif? = 

<5i< = TVfC (3) 



determined by a single grassmann parameter ^ corresponding to a single 
supercharge. We know that the continuum N = 2 Wess-Zumino model pos- 
sesses four such supercharges corresponding to two independent two com- 
ponent Majorana charges. Thus we might expect that the lattice model will 
admit three further transformations which become invariances as a — > 0. 
The complex form of the bosonic action immediately suggests three further 
bosonic actions which will differ from each other by terms which become 
total derivatives in the continuum limit. These are 



^(2) 


= D,4 


> - W'{(t)) 




= D,4 


> - W'(^) 


^(4) 


= 


> + W'(^) 



(4) 

Let N"^ be the (real) two component vector corresponding to the complex 
field ri^'^\ Under differentiation it generates a new fermion matrix 



ap _ 



Using the arguments of the previous section we can now write down a new 
lattice action S^^^. 

5(2) = lr,(2)r^(2)+;^M(2)a; 

where x ^''^d u are new anticommuting spinor fields. S^^^ will, of course, 
possess a new supersymmetry invariance involving now not the vector N 
but N. Furthermore it is easy to see that M*^^^ = ij^Mij^. Hence these 
two lattice theories generate (up to total derivative-like terms) the same 
continuum action. Indeed, if we make the identifications 

tp = i^sLO 

ijj = i-fsx (5) 
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we can see that the original lattice action has a second approximate super- 
symmetry given by 

62^? = 

62i^t = ilfN^ (6) 

The variation of the action under this second supersymmetry involves the 
supersymmetry variation of terms which vanish as total derivatives in the 
continuum limit. On the lattice these terms will be of order ga^ with g a 
typical interaction coupling. Hence, at least in perturbation theory such a 
term would constitute an irrelevant operator and the continuum limit should 
exhibit this second supersymmetry. One might worry that the presence of 
such a SUSY- violating term in the bare lattice action might lead to relevant 
breaking terms in the long distance effective action. However, it is not 
possible to write down any such counterterms which simultaneously preserve 
the one exact SUSY. Thus, the existence of a subset of the full SUSY in the 
lattice model is indeed sufficient to protect the broken supersymmetries so 
that no fine tuning is required to achieve the full symmetry in the continuum 
limit. 

Turning to r/^^^ we can see that it generates yet another fermion matrix 
of the form 

= 5 

where the vector Q° again carries the real and imaginary parts of ri^^\ 
Again, M^^) may be expressed in terms of the original M 

which proves that an action based around rj^^^ will once again constitute a 
lattice theory of the continuum Wess-Zumino model with yet another super- 
symmetry. In terms of the original fermion fields this third transformation 
will yield another approximate invariance of the original action 

Ss^t = 

Ssi^f = Qti (7) 

The final approximate invariance can be derived similarly from tj*^^-* (or its 
real vector form ) and yields the transformations 

5^x1 = iftii 
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5,i:t = ilfQk (8) 

Thus far we have again assumed that the variation of the fermion matrix 
under these supersymmetry transformations is zero. However, the simple 
proof we gave in the previous section for the absence of such a term in 6S 
does not hold when the variation of the field x involves non-trivial gamma 
matrices acting on ^l; or ip. If we examine the general structure of such a 
variation we find that it has the form (we suppress spacetime indices which 
play no essential role) 

dxl^dx'y 

where 0, 9 represent either ^ or ip. This can be seen to be the trace of a 
product of a symmetric matrix (the term involving derivatives of /) with the 
gamma matrix F and the antisymmetric matrix formed by the product of 
the 9 terms. Thus, for F = 71 or F = 72 this is the trace of an antisymmetric 
matrix and is hence zero. For F = ^73 the resultant matrix is now symmetric 
but the trace can be shown to still vanish as a consequence of the Cauchy- 
Riemann conditions applying to the derivatives of /. 



4 Ward Identities 

4.1 Quantum Mechanics 

The invariance of the quantum mechanical lattice action under the discrete 
supersymmetry transformation eqn. |2| leads to a set of Ward identities con- 
necting bosonic and fermionic correlation functions. We can derive these 
following the usual procedure by adding a set of source terms to the action 
and carrying out an (infinitesimal) supersymmetry variation of the fields. 
Since the partition function, measure and action are all invariant under this 
change of variables we immediately derive the result 

6Z = = J Di^Di^Dxe-^^-^-''^^^ ( J.<5ix + 9.5i'^^ 

Indeed any derivative of this expression with respect to the source terms 
(which are set to zero at the end) is also vanishing. Thus we are led imme- 
diately to the first non-trivial supersymmetric Ward identity 

l^iij,) + {NiXj) = 
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relating the fermion correlation function to one depending only on bosonic 
fields. Notice also that in the continuum limit there will be a second set of 
Ward identities following from the second invariance given by the variation 
S2. _ _ 

To perform a simulation of this model we will replace the integral over an- 
ticommuting fields ip, ip hy one over a (real) pseudofermion field x whose ac- 
tion SpF = x^iM^M)-^x yields the same fermion determinant det(M(x)). 
Consider now the generalized partition function Z{a) where 

Z{a) = J DxDxe-'''^^'''^^ (9) 

This allows us to write down a simple expression for the mean action in- 
cluding the pseudofermions 

^ _ dlnZ{a) 



da 

We will from now on restrict ourselves to lattice actions which derive from 
a field Ni of the form 

Ni = DijXj + MijXj + gxf 
In this case a simple scaling argument allows us to rewrite eqn. ^ as 
Z{a,9)=a-''/^Z{l,g') 

where g' / g = a' 2 and is just the total number of degrees of freedom we 
integrate over. Hence we find the following expression for the expectation 
value of the total action including the pseudofermions 

N 1-Q d 

The second term on the right vanishes by virtues of the fact that the partition 
function does not depend on gf - as guaranteed by the existence of the Nicolai 
map. Thus we see that the mean action (with a = 1) merely counts the 
number of degrees of freedom including the pseudofermions. Furthermore, 
since the existence of the Nicolai map implies a supersymmetry we can also 
regard the value of the mean action computed in the simulation as an order 
parameter for supersymmetry breaking - if we find it depends on coupling 
and differs from its value for the free theory we know that supersymmetry 
has been broken. 
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4.2 Wess-Zumino Model 

The analysis of the previous section carries over to the Wess-Zumino model 
with the appropriate interpretation of the index and field content. Thus 
we expect the mean lattice action to be equal to the number of degrees 
of freedom < S >= 2L^ for a lattice of linear size L (the two counts the 
two real degrees of freedom at each lattice point in either boson or fermion 
sector). Likewise we expect the Ward identity based on the variation 5i to 
be exact for arbitrary lattice spacing. 



Similarly we expect the 
a ^ 0. 

= 
= 
= 



(CV'f ) + (iv/x-) = 



(10) 



following three Ward identities to be satisfied as 

(irr^^^) + (^I's'Q]^) (11) 



5 Numerical Results 

To check these conclusions we have chosen to simulate the model for W'{(j)) = 
nKp+gcf^. We have used a hybrid monte carlo algorithm Q to handle the in- 
tegration over the pseudofermion fields. In order to reduce the computation 
time for large lattices we have implemented a refinement of this algorithm 
using Fourier acceleration techniques. Details are given in |p and more re- 
cently [|l^]. In the latter paper we show that the autocorrelation time for 
SUSY quantum mechanics is drastically reduced - the dynamical critical ex- 
ponent z is reduced from z ~ 2 for the usual HMC algorithm to z ~ with 
fourier acceleration. In the Wess Zumino case the gains are also large. 

Weak Coupling 

In order to compare our results with other continuum and perturbative 
calculations we simulated the model initially at zero and small coupling g. 
We show data for m = W g = and g = 3 obtained from 1 x 10^ HMC 
trajectories at L = 4, and L = 8, 2 x 10^ HMC trajectories at L = 16 and 
2 X 10^ HMC trajectories at L = 32. To take the continuum limit we imagine 
holding the physical size of the lattice fixed at unity (we are neglecting 
finite size effects sinhce our bare masses are relatively large). This allows 



12 



L 


g = 0.0 


g = 3.0 


4 


32.01(4) 


31.93(6) 


8 


127.98(6) 


127.97(7) 


16 


512.5(3) 


512.0(3) 


32 


2048(1) 


2046(3) 



Table 1: Mean total action vs lattice size 



us to extract the lattice spacing a = Since our lattice action contains 
only dimensionless quantities the bare physical couplings g and m must be 
translated to bare lattice quantities = g/L, = m/L in the lattice 
action. The continuum limit is then reached by simply taking L — ^ oo. 

Table 1. shows the mean action as a function of lattice size for both 
g = and g = 3. As is evident the mean action is close to the predicted 
value of 2L^ consistent with a non-breaking of SUSY (this is expected since 
the Witten index for this model A = 2). 



0.02 




5 10 15 

t 



Figure 1: Boson Correlator at L = 16 and m = 10.0, g = 3.0 
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To extract information on the spectrum of the model we have studied 
zero momentum correlation functions which are given by averaging the fields 
transverse to the direction of propagation. 

hi' 

and 

On account of the periodic boundary conditions we expect the boson corre- 
lator G^{t) =A{t-L/2) where A is a symmetric function of its argument. 
Conversely the fermionic correlator can be expected to take the form 

G^p{t) = k{I^pB{t - L/2) + \^ipC{t - L/2)) 

where B{x) and C{x) are symmetric and antisymmetric functions of their 
arguments, A is a numerical coefficient and 7* is the gamma matrix appropri- 
ate to the t-direction. For large x we expect a single mass state to dominate 
in which case A{x),B{x) — > cosh (m^x) and C{x) — > sinh(m^x). These 
latter functional forms were found to yield good fits over the whole range of 
parameters studied. The parameter nig corresponds to the massgap of the 
model expressed in lattice units. To convert this value to physical units we 

merely have to divide by the lattice spacing a, nig = . 

Fig. I and fig. | show Gf^(t) and G^p{t) for L = 16, 5 = 3 a = 
(3 = 2 and t lying along the 1-direction. This choice of time direction 
implies that the fermion correlator will be purely diagonal with Gii{t) = 
cosh (m^(t - L/2))-hAsinh (m^(t -L/2)) and GfaW = cosh {ni^{t -L/2))- 
Asinh(m^(t — L/2)). For weak coupling we find that the numerical value 
of A extracted from the fit is consistent with unity which would be expected 
for a free theory as a — > 0. Notice that although the lattice action does not 
satisfy refiection positivity there is no sign of a problem in the correlation 
functions at weak coupling. 

In Table 2. we show the results for the massgaps in physical units rUg 
as a function of the lattice spacing. It is clear that the boson and fermion 
masses are degenerate within statistical errors and increase smoothly with 
decreasing lattice spacing. 

Fig. m is a plot of physical (fermion) mass nig{a) extracted from the 
simulations as a function of lattice spacing. For small g/ni = /m^ we 
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L 


rriB 


mp 


4 


5.09(2) 


4.95(8) 


8 


6.52(2) 


6.44(5) 


16 


7.76(4) 


7.75(6) 


32 


8.29(19) 


8.33(30) 



Table 2: Physical massgaps nig vs lattice size for m = 10.0 and g = 3.0 



expect perturbation theory to provide a good approximation. The one loop 
result for the massgap is 



which yields m^^'^^ = 9.65 for g = 3.0. Notice that since this theory is finite 
there is no need to introduce a scale dependent renormalized mass - the 
physical massgap of the theory is a finite function of the bare parameters in 
physical units. It is encouraging that reasonable extrapolations of mg{a) to 
a = are consistent with the one loop result. These numerical results are 
also consistent with ones which were previously obtained using a stochastic 
approach based on the Nicolai map |11|. 



To understand whether the continuum limit will describe an = 2 
supersymmetry we have also checked the four Ward identities written down 
in the last section. We again choose to average the correlations transverse 
to a chosen t-direction (t = 1 as before). Each Ward identity then yields two 
independent relations between components of boson and fermion correlators. 
Fig. § shows a plot of —G22{t) and < X2{^)N2{t) > versus time t. The first 
Ward identity requires the sum of these two curves to vanish - clearly to a 
very good statistical accuracy the numerical data support this conclusion. 
The first Ward identity also predicts a relationship between Gf|(t) and < 
xi{Q)Ni{t) > which we also observe to be true within (small) statistical 
errors. Thus, as expected, the existence of the exact SUSY eqn. |^ leads 
to a Ward identity relating correlation functions which we observe to be 
accurately satisfied on the lattice. 

We have also examined the other Ward identities corresponding to the 
other three broken symmetries - fig. ^ plots G22{t) vs < xi(0)A^i(t) >. 
Again, if the 2nd Ward identity were to hold exactly the sum of these two 
curves would again vanish ~ and it appears that the data are consistent with 
this. Indeed, we have found that each of these three Ward identities is also 
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satisfied within statistical error at this (weak) coupling. The explanation 
for this seems to lie in the magnitude of the symmetry breaking - as we have 
argued the breaking effects come in at 0{ga^) = g/L^ yielding corrections 
to the broken Ward identities which are too small to be resolved over our 
statistical errors. 

Strong Coupling 

Having checked by explicit simulation that this lattice model appears to 
possess the correct supersymmetric structure at weak coupling we have ex- 
tended our simulations to strong coupling. This allows us to probe directly 
the non-perturbative structure of the theory. Classically the model has two 
vacua - corresponding to </) = and 4> = —f^/ 9- These vacua are separated 
classically by a barrier of height {m/Ag)^ . For small g we are effectively 
confined to the i;^ = well but for large g/m we would expect both vacua 
to be sampled. In addition since the terms in the action that violate re- 
flection positivity are proportional to g we might wonder whether a sensible 
continuum limit even exists for strong coupling. 

We have examined this issue by simulations at m = 5 and g = 2.5, 
g = 5.0 and g = 10.0 for lattices from L = 8 through L = 32 as before 
and with similar statistics. The choice of a smaller bare mass parameter 
m reduces the barrier height and allows our simulation to more effectively 
tunnel between the two classical vacua. For large g we were forced to refine 
our Hybrid Monte Carlo scheme to eliminate problems stemming from large 
pseudofermion forces occasionally encountered in the vicinity of such tun- 
neling configurations of the boson field. Essentially an entire trajectory is 
abandoned as soon as a force component larger than some threshold is seen 
- the trajectory is restarted with new momenta. This process increases the 
autocorrelation time of the algorithm but for the parameters at which we 
performed simulations the effect was not overly severe. 

The results for g = 2.5 are similar to those obtained in the previous 
section and will not be examined further. In all cases we have observed 
that the mean action < S* >= 2L^ independent of g and m. This is again 
evidence that supersymmetry is not broken even outside of perturbation 
theory. We have also seen that typical configurations extend over a region 
of field space encompassing both classical minima - indeed we have found 
that < Re0 >= —m/2g very accurately. 

The correlation functions G'^(t) and G22{t) for g = 5.0 [glm = 1) are 
shown in fig. ^ and fig. 0. We again choose time along the 1-direction. 
The boson is again accurately fitted by a simple hyperbolic cosine function 
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L 


rriB 


mp 


8 






16 


4.35(7) 


5.5(4) 


32 


6.0(2) 


4.4(7) 



Table 3: Physical massgaps vs lattice size for m = 5.0 and g = 5.0 



and yields a physical mass of = 4.35(7) at this lattice spacing. The 
fermionic correlator is a little more complicated - at this coupling the fits 
favor a signal which is predominantly given by a hyperbolic sine function 
with a small admixture of hyperbolic cosine (typically A ~ 4/5). The use 
of a three parameter fit yields a larger error in the fermion mass estimate - 
mp = 5.5(4). The results for all the lattice sizes are summarized in table. 
3. 

It is not clear whether the discrepancy between boson and fermion mass- 
gaps is significant or merely reflects the large errors in the fermion mass 
determination. More interestingly, the gaps in the table for L = 8 arise 
because it was not possible to extract a mass from the small lattice L = 8 
- the signal descends into noise after just one timeslice. This is not true 
for smaller g and may indicate a problem with reflection positivity at this 
lattice spacing. A similar problem occurs for L = 16 when g is increased to 
g = 10.0. Fig 1^ shows a plot of the bosonic correlator there. We conjecture 
the oscillations visible in the signal are a signal for a mass spectrum which 
is not real positive. This might indicate that the problem could indeed be 
attributed to the lack of reflection positivity in the lattice action. However, 
even if this were the case, the problem appears to diminish with lattice spac- 
ing - the correlators for L = 32 at this same coupling g = 10 exhibited none 
of these problems and allowed fits for both boson and fermion massgaps 
niB = 4.7(1) and mp = 4.9(7) for L = 32. Interestingly, in the region of pa- 
rameter space where the two-point functions show this oscillatory behavior 
we have also observed that the sign of the fermion determinant may fluctu- 
ate also. In practice the sign changes were relatively infrequent and their 
effects could be taken into account by reweighting the measured observables 
in the usual manner. However, this effect also disappeared with decreasing 
lattice spacing at fixed coupling. 

Similar results were obtained at other values of the bare parameters. 
Hence we speculate that, while problems associated with a lack of reflection 
positivity may be evident on coarse lattices, these effects disappear with 



17 



decreasing lattice spacing. Thus a well-defined continuum limit may be 
defined for all finite g. 

Finally we have examined the Ward identities. Again, the presence of 
an exact symmetry yields a relationship between boson and fermion cor- 
relators for arbitrarily large g which is exhibited in fig. ^ which shows the 
same correlators deriving from the first Ward identity now for m = 5, g = 5 
and L = 16. The middle curve (diamonds) shows the sum of the two con- 
tributions which is seen to be consistent with zero for all t within errors. 
Contrast this with fig. |l^ which exhibits the bosonic < xi{0)Qi{t) > and 
fermionic Gfi {t) contributions to the third Ward identity. The middle (dia- 
mond) curve is no longer zero and indeed shows a marked variation with t. 
Similar effects are seen in the fourth Ward identity. 

6 Discussion and Conclusions 

In this paper we have studied a lattice version of the two-dimensional Wess- 
Zumino model with N = 2 supersymmetry. The lattice action we use was 
first derived in , Q and follows from a discretization of the continuum Nico- 
lai map for the model. We have rederived it in a different way by requiring 
that the lattice field theory model exhibit a single parameter SUSY-like in- 
variance. This approach has the advantage that is allows us to identify the 
other broken invariances which would yield a full N = 2 SUSY in the naive 
continuum limit. From the form of those transformations we have derived a 
set of Ward identities which would be satisfied in the continuum limit. We 
furthermore argue that the presence of one exact symmetry (together with 
the finiteness of the continuum theory) guarantees that the full symmetry 
is restored without fine tuning in the continuum limit. 

These conclusions have been checked by an explicit numerical simulation 
of the Euclidean lattice theory in which the boson and fermion massgaps 
and a set of supersymmetric Ward identities were computed at a variety 
of lattice spacings. We utilized a Fourier accelerated Hybrid Monte Carlo 
algorithm to handle the fermionic integrations. 

At weak coupling we were able to extract boson and fermion masses 
and verify their equality within statistical errors. We also found that all 
the Ward identities were satisfied to high precision. We have argued that 
the small magnitude 0{ga'^) of the symmetry breaking effects places the 
corrections within the statistical noise inherent in our calculation. Most 
importantly, the numerical results show no sign of any problems stemming 
from the lack of reflection positivity in the lattice action. 
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At strong coupling we found difficulties in extracting masses and inter- 
preting the theory for coarse lattices but, at least for the parameters we 
studied, these effects seemed to go away on finer lattices. Our simulations, 
while efficiently sampling the classical vacua of the model, show no evi- 
dence for supersymmetry breaking - the mean action remained at 2L^ and 
the Ward identities corresponding to the exact symmetry were still satisfied 
within errors. However at strong coupling we did observe clear corrections 
to some of the other approximate Ward identities. 

Of course the interesting question is whether one can generalize any of 
these ideas to gauge models in higher dimensions. A Nicolai map is known 
for the continuum N = 2 super Yang Mills model in two dimensions [^] 
(indeed it can be obtained by dimensional reduction of a map for = 1 
super Yang Mills in four dimensions). Unfortunately, a naive transcription 
to the lattice is problematic since the map utilizes an explicit noncompact 
formulation for the gauge field. Replacing continuum derivatives by finite 
differences as for the Wess Zumino model would then lead to an action which 
was not gauge invariant. 
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Figure 5: Fermionic and Bosonic Contributions to 2nd Ward identity for 
m = 10.0, g = 3.0 
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Figure 6: Bosonic Correlator m = 5, g = 5, L = IQ 
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